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Abstract 

We consider stochastic inflation in an interacting scalar field in spatially homogeneous accelerating 
space-times with a constant principal slow roll parameter e. We show that, if the scalar potential 
is scale invariant (which is the case when scalar contains quartic self-interaction and couples non- 
minimally to gravity), the late-time solution on accelerating FLRW spaces can be described by a 
probability distribution function (PDF) p which is a function of p/H only, where p = p{x) is the 
scalar field and H = H{t) denotes the Hubble parameter. We give explicit late-time solutions for 
P Pooip/H), and thereby find the order e corrections to the Starobinsky-Yokoyama result. This 
PDF can then be used to calculate e.g. various n—point functions of the (self-interacting) scalar field, 
which are valid at late times in arbitrary accelerating space-times with e = constant. 
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I. INTRODUCTION 


Q 


In 1986 Alexei Starobinsky [Ij realised that a classical stochastic theory can be used to 
calculate infrared correlators of interacting scalar helds on De Sitter space. The stochastic 
formalism was then used to calculate the late-time state of a self-interacting scalar held j^. 
Almost 20 years later Tsamis and Woodard proved that stochastic formalism correctly cap¬ 
tures the leading logarithms on De Sitter background to all orders in perturbation theory, and 
therefore can be used to efficiently resum the perturbation theory on De Sitter in scalar inter¬ 
acting theories. Miao and Woodard jd] and Prokopec, Tsamis and Woodard have shown how 
to generalize the results of to include light fermionic and (massless) Abelian gauge helds, 
respectively. Since massless fermions and Abelian gauge helds couple conformally to gravity, 
they cannot be treated as classical. The authors of 4] and realized that these passive helds 
(light fermions and Abelian gauge helds) can be integrated out to obtain an ehective interacting 
scalar theory, which can be then stochasticized Iw the usual Starobinsky’s procedure. Explicit 
two-loop calculations in scalar electrodynamics [6 71 vindicate the stochastic inhationary pre¬ 


dictions of . Despite of some notable attempts |3|, [8 


10| . it is fair to say that at this moment 


it is not known how to rigorously stochastize (perturbative) quantum gravity on De Sitter or 
on inhationary space-times. 

Recently Gihyuk Cho, Cook Hyun Kim and Hiroyuki Kitamoto [ffi, 121 showed that the 
correct way to stochasticize light interacting scalar helds on accelerating spaces (dehned by 
0 < e = —H/H^ < 1, where H is the Hubble parameter), introduces an important order 
e-correction to the friction in the naive stochastic equations. In this letter we make use of 
Cho-Kim-Kitamoto’s result to solve for the late-time solution of the held. 

We show that scale invariant interacting scalar theories admit at late times a solution in 
which the PDF p = p(^/H) = exp[—Ves{^/H)], where Ves{^/H) generalises the Starobinsky- 
Yokoyama solution 2|. Recall that the Starobinsky-Yokoyama solution is of the form, 
PsY = Poexp(—n), with v = [8n‘^/ {3 H‘^)]V{ip), where V{ip) = A(p'^/4! is the tree-level po¬ 
tential. Strictly speaking this solution applies to De Sitter space only. In this letter we show 
that in accelerating spaces with a constant e parameter, the late-time solution of Starobinsky 
and Yokoyama psY gets modihed. Furthermore, we explicitly calculate the e-induced corrections 
to the PDF. Next, we show that only some classes of potentials - namely those that exhibit 
scale invariance - permit simple late-time solution in the form p —>■ poa{p/H). Curiously, these 
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are precisely the potentials that coincide with those studied in Ref. [13| that admit classical 
solutions of the form, 0(f) = 'ijjoH{t), where i/jq is some (calculable) constant. In Ref. 13| the 


one-loop-effective potential was calculated for a background held of the form, 0(f) = 'ijjoH{t). 
The conclusion of Ref. was that two types of quantum effects break the classical scale 
invariance: (1) the effects induced by the counterterms in the process of renormalization and 
(2) the effects resulting from regulating the infrared sector of the scalar held (whose naive 
Bunch-Davies vacuum contains power low divergences in the infrared). Remarkably, stochastic 
analysis presented in this letter shows that, when perturbative contributions are resummed, 
the scaling symmetry of the classical theory gets restored, constituting the main result of this 
letter. 


II. THE MODEL AND ITS STOCHASTIC ANALYSIS 


Consider a scalar held 0 whose dynamics is implied by the (classical) action, 
S'[0] = f 




( 1 ) 


where is the metric tensor, g^^’^ its inverse, the metric signature is sign[ 5 f^,^] = (—1,1,1,1), 
R = R{g^v) is the Ricci scalar, and R(0) is the (tree level) scalar potential. When gravity is 
taken to be non-dynamical and when the scalar held is assumed to be light, | [cf^R/cf0^]<^=o| "C 
|.^| -C 1/12, the scalar theory ([1]) on De Sitter space and in the infrared (on super-Hubble 
scales) is well approximated by Starobinsky’s stochastic inhation [l, 2|. In this letter we study 
the late-time dynamics of the scalar held by making use of the modihed stochastic inhation hrst 
obtained in jlll. Il2|. To make it suitable for stochastic treatment, the coupling constants f and 
the corresponding couplings in the potential V (0) are the renormalized couplings evaluated at 
the physical momentum scale /2 ~ fihfo; where ju 1 and ffg is the Hubble parameter at a 
suitably chosen initial (sufficiently early) time. 

Here we work in a hat FLRW space-time, for which the metric tensor can be written as. 


= —dt^ + of{t)dx- dx , 


( 2 ) 


where a = a{t) is the scale factor. In terms of the scale factor, one can dehne the Hubble 
parameter, H{t) = a/a and the principal slow roll parameter, e = —H/H’^. The Ricci scalar is 
then, R = 6(2 — e)H‘^. In De Sitter space, if = 0, such that e = 0 and R = 12H^ is constant. 
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A. Stochastic theory on De Sitter space 


Stochastic inflation on De Sitter can be obtained by the following simple procedure (which 
is not a derivation; for a derivation see js, 5|). The action ([T]) implies the following operator 
equation of motion for 0 , 


(j)(t, x) + 3H(f) - -(j) + D'(0) + ^R(j) = 0 , 


(3) 


where V = dV/d(j). Here we are interested in the dynamics of the infrared part <h(t, T) of the 
held <p. In the limit when the scalar held is light, |m^| = = 0) + the 

scalar held dynamics is correctly described by the following simple stochastic equation of the 
Langevin type, 


$ + 


!/'($) - 



(4) 


where the couplings in V (<h) and ^ are to be evaluated at the ultraviolet cutoff scale of the 
theory, given hy fi = fiHo (p 1 ), and $ 0 (^ 5 is the solution of the homogeneous linear 
equation, 

<f)o(t, T) + + [6a2-e)H^ + H"(0)] $0 = 0 . (5) 

CL 

This linear equation can be solved by decomposing the held in terms of the mode functions 
u{t, k) and u*{t, k) as. 


$o(t,T) = 


d^k 


Q{fiHa — k) u{t, k)d{k)R^'^ + u*{t, k)d'{k)e 


-ik-x 


( 6 ) 


where d{k) and d){k) are the annihilation and creation operators and /i -C 1 is the ultraviolet 
cutoff scale. The operator d{k) annihilates the free vacuum state |D) of the theory, d{k)\Q) = 0, 
and cd{k) creates a particle of momentum k out of the vacuum. Applying successively a){k) on 
|f2) generates the Fock space of the theory. 

On De Sitter space and in the massless limit (e = 0, = 0), the properly normalized 

mode functions (which obey the canonical Wronskian, W[u,u*] = i/o?)-, are given by. 


u{t, k) = 


H 


ik 




exp 




u*{t, k) = 


H 


ik 


vW L 


exp ( — i 


aHJ 


(7) 


In the more complicated case, when the scalar is massive, the mode functions can be expressed 
in terms of the Hankel functions of the first and second kind. From the solutions ([7]) it is 
clear that the first term in the square brackets dominate in the infrared (on super-Hubble 
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scales), when k/a H, such that the mode functions in $0 in Eq. (|6]) can be approximated 
by, u{t, k) = k) ~ H/\/2k^. With this approximation in mind, one can show that the held 
operator $0 becomes a classical stochastic held, in the sense that its canonical commutator 
vanishes, [<ho(t, x), $o(t, ^)] = 0, and its correlator is Markovian, 




(pa)^(/ii/a) sm[fj,Har] 


187r2 


^Har 


\u(t, k = fiHa)\‘^S(t — t'), r = ||x—^11, (8) 


where ^{t,x) = <ho/(3-ff). The formula (jH]) simplihes further when one takes account of the 
infrared approximation for the mode functions, \u{t,k)\‘^ —)■ l/[2H{fia)^] and in the limit of 
spatial coincidence, when sm[^Har]/{^Har) 1. In this case the noise correlator ([H]) simplihes 
to, ~ [H/{367i^)]6{t -t’). 


B. Stochastic theory on accelerating spaces 


Cho, Kim and Kitamoto recently showed [nl,Q that, in order to obtain the stochastic 
theory applicable to accelerating spaces with constant e (for which 0 < e < 1), it is convenient 
to rewrite the operator equation ([3]) as 


dj^(j){t,x) + (3 - e)c?Ar0 


{aHf 


+ 




= 0 , 


(9) 


where N = ln(a) is the number of e-foldings. The corresponding stochastic theory is then 
obtained by dropping the second derivative terms in ([9]) and adding a stochastic force ^ as, ^ 


aAr$(iV,f) + 


!/'($) - l/"(0)<h 
(3 - e)i/2 


$0 

(3 - e)H 




where the stochastic force ^ obeys. 




(/ia)^i/^(l — e) sin[nHar] 
(3 — e)227r2 nHar 


u{t, k)\'^6{t — t '), 


( 10 ) 


( 11 ) 


^ There is no unique way to split the potential into the one that gives rise to the stochastic force and the 
interactions. In Eq. we use the standard procedure, and include the quadratic terms in the action into 
the equation of motion for the free field and all cubic and higher order terms into the interactions. This 
split is supported by the well known fact that large masses suppress stochastic force, hence it is natural to 
include them into $o- In the Appendix we consider another possibility, according to which mass terms are 
split between the interactions and the free field. That choice is motivated by the desire to retain the cutoff 
independence of the De Sitter case. 
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and u is the mode function satisfying, 


dl + e- — + a^ [i/"(0) + 6^(2 - e)H^] ) [au{7], k)] = 0. 

Qj / 


( 12 ) 


For a constant epsilon, a"/a = (2 —e)/[ 77^(1 — e)^] and assuming m^/H = const, (we will see 
below in which situations this approximation holds), we get 

^2 (2-6)(1-60 1 




ii-ey 


+ av ''{^)) [au{ri^ k)] = 0. 


(13) 


The Bunch-Davies (or Chernikov-Tagirov) vacuum solution ^ for the mode functions is (see e.g. 


Refs. 


131 and 


Q). 


u{7],k) = - 

a 

u*{v,k) = - 


TT 


4(l-e)-H 




k 


vr 




4(l-e)-H 


k 


^ = 7 + 
4 


1 , (2-6)(1-60 R"(0) 


(1 




(14) 


where "H = d\ln{a{r]))]/dT] = aH is the conformal Hubble rate. These solutions are correct 
provided V "/are time independent (in slow roll inflation one allows for a weak dependence 
of p on time, in the sense that \i'/v\ if, which amounts to an adiabatic approximation). 
Since mode functions in (ITTll are truncated at k/a = /iff <C ff, one can approximate the 
Hankel functions in flTTll with their small argument expansion. ^ When the leading term in this 
expansion is inserted into (ITT]) one obtains, 

h) sinl/iffarl ^ 


where 


A = 


/i 


(3-6)2 
3-2-(l_g) 


fiHar 
2^22'^-3r(z/)2 


vr-" 


(16) 


such that, when 1 / ^ 3/2, stochastic force becomes dependent on the cutoff scale /i. If we are 
interested in the cases when |e| -C 1, 12|,^| -C 1 and |R"(0)| ff^, then v 3/2 (|i/ —3/2| -C 1), 
the dependence on the UV cutoff is weak, and stochastic treatment applies. One could argue 
that the better choice for the UV cutoff in Eq. (jSj) is fuv = f(1 ~ e)aH. Indeed, for this choice 
the argument of the Hankel function becomes /r, which equals to 1 when /r = 1, signifying the 


^ The mode functions for the most general homogeneous pure Gaussian state can be written as, a{k)u{g, k) + 
(3{k)u*{g, k), with |ap — |/3p = 1- 

3 The relevant small argument expansion of the Hankel functions is, H^\z) = [H^\z)\* ~ \r[vY/' k\{ 2/zY'' + 

C>(z-2'^+2,2°, (3 [z 2] = 0, > 0, 5[d = 0). 
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scale at which mode functions start coupling strongly to the background, and the process of 
large (non-adiabatic) particle creation sets in. For that choice of the cutoff, A in Eq. flT6|l would 
rescale as, A A{1 — , retaining thus the same dependence on the cutoff scale p, thus 

not changing flTBD in any significant way. 

Just like in De Sitter space, the stochastic theory ffTUHTT]) can be recast as a Fokker-Planck 
equation ^ for the probability distribution function (PDF) p = p{N{t), ip{x)) (for the stochastic 
field (p{x) which represents the quantum field <F(f,x)), 


dNp{N, (p(f)) = ^^dlp + ’ (1'^) 

where we made use of dt = Hd^, which follows from N = ln(a). Starobinsky and Yokoyama 
told us how to solve the Fokker-Planck equation in the case of a self-interacting massless scalar 
field on De Sitter space, for which D(0) = [A/4!]0^. Their solution can be easily extended to 
the case of a light scalar, even though in that limit the solution picks up a slight dependence 
on the cutoff scale p. 


C. Stochastic late-time solution 


In order to study asymptotically late-time solutions of flTTl) in an accelerating space-time 
with constant epsilon, it is convenient to transform from {N, (p) to the variables {N, (p = p}/H), 
in which case, 

—)■ -\- e(pd(p . ^ H ^d(p . (18) 

Furthermore, from j dipp{N,ip) = f d(pp{N,(p), it follows that p = Hp. ^ With these in mind, 
Eq. (ITTll becomes, 

dNp{N, Cp) + ed^ [pp{N, Cp)] = + ^5^[pY'(<^)], (19) 


where V {x) is defined by. 






( 20 ) 


A closer inspection of flT^ implies that, if we are to have a time independent solution at 
late times, then time dependence must disappear from the right hand side of (IT^ . i.e. V{(p) 


I thank Hiroyuki Kitamoto for pointing out to me that this rescaling of p significantly simplifies the late-time 
solution for the field PDF. 
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must be time-independent. This is the case when the scalar held is invariant under scale 
transformations, ^ i.e. when V{(j)) = [A/4!]0^. In this case, 


r(v) = 


( 21 ) 


For example, neither cubic interactions, Asip^/S!, nor tree mass terms, m‘^ip'^/2, are not allowed 
in the potential, because it would generate the terms, [m//2 + {X^/QH)(p^, in V, which 
depend on time through 

With these comments in mind, the late-time limit of Eq. fll9l) becomes. 


= jdip^ + 


( 22 ) 


This equation governs the late-time PDF of the scalar held in stochastic inhation. Next, it is 
convenient to rewrite Eq. (|2^ in the form, 

2 


dpppoo T d(p ^ 1 ^eff 


(3-e)2l L 


f>/~\ (3-^)^.2 


(23) 


Then the following transformation, = po exp[—n(<p)] yields the diherential equation for n(<p). 


+ [-' - vUV = 0 . 


(24) 


The obvious (and correct) solution is. 


v{if) = lieu)*;), 


(25) 


where VeE{p>) is given in Eq. (|23|) . In conclusion, the evolution of an interacting scalar held in 
accelerating spaces gets modihed such that - due to the deviation from De Sitter exponential 
expansion - at late times the held acquires a negative, time dependent, mass-squared, 

= -(3-e)eif2 (26) 

Since our held already has a (time dependent) mass-squared, R = 6^{2 — e)H‘^, the total mass 
becomes, 

[6^(2-e)-(3-e)e]i/". (27) 

® By the invariance under scale transformations we mean the invariance of the classical action under the 
transformations of the type, g^ti, —> ^ where P is an arbitrary constant. It then follows that 

the Ricci scalar transforms as, R —>■ P“^i?, \J—g — >■ and the action ([T]) is invariant if F —P^'^F, 

which is the case when F = Ay'^/4!. 










Depending on the sign and magnitude of this mass-squared can be either positive or nega¬ 
tive, indicating that deviations from De Sitter can be responsible for generation of symmetry 


breaking. A complete understanding of this interesting ques 


discussion of symmetry restoration on De Sitter see Refs. 


15 


ion requires further study (for a 


131 ) 


The PDF p{lp) that we obtained in this letter can be used to calculate asymptotically late¬ 
time coincident held correlators in accelerating spaces. In order to obtain the corresponding 
renormalized correlators, the following formula can be used. 


{(p {x)) ~ ((^ )stoch = / = hr / dipp{ip)(f 


-2n 


(28) 


whereby p is assumed to be suitably normalized, j d<pp{(p) = 1. The relevant integrals can 
be expressed in terms of the conhuent hypergeometric functions and the modihed Bessel func¬ 


tions 


id- 


III. CONCLUSION AND DISCUSSION 


In this letter we have constructed the late-time solution for an interacting scalar held the¬ 
ory on accelerating spaces with a constant principal slow roll parameter, e = The 

most general tree level potential for which one can construct an exact solution contains quartic 
self-coupling and a non-minimal coupling to gravity. Our solution fl2^ shows that, the ehec- 
tive potential that governs the PDF at late times acquires a negative, order e mass-squared, 
Am^ = — (3—e)eih^, see Eq. (|26|) . thus ehectively increasing the size of held huctuations. If the 
total mass-squared becomes negative, that would suggest a symmetry breaking induced by the 
expansion. It would be of interest to understand the dynamics of such a symmetry breaking. 
Even in the simpler De Sitter background the problem of symmetry breaking/restoration has 


not fully been understood 


15 


17l |. The reason why we were able to hnd an exact late-time solu¬ 


tion for the held PDF is because the tree-level potential is invariant under scale transformations. 
The resulting ehective potential (|25il . (|23|l inherits the scale invariance of the tree-level. 

An interesting question is how our results get modihed when one includes into the tree 
potential terms that break scale invariance, such as a mass term or a cubic self-interaction. 
It is easily seen that in that case the Fokker-Planck equation flTTl) does not admit late-time 
scale invariant solutions simply because in this case the rescaled potential V contains time 
dependent terms V D (f)^] cx A 3 (p^/l 6 II(t)] oc It would be of interest to 
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construct adiabatic solutions in that more complicated case. We suspect that is possible only 


I2I. 


when e is small, i.e. e <C 1. For a discussion of that question see Ref. 

One annoying feature of our solution (I2B]) . (125]) is the cutoff dependence of our late-time 


solution for the PDF p{ip) , which is hidden in the /r-dependence of the stochastic force param¬ 
eter (1161) . A oc That dependence is weak when z/ ~ 3/2, which is the case during slow roll 

inflation and when the held mass is negligibly small. There is a formal way of getting rid of this 
fi dependence, which is discussed in the Appendix. The idea is to split the mass term into two 
parts, such that the stochastic force is equal to that of the massless held in De Sitter limit. Just 
like in the De Sitter case, the dependence on p then disappears. The resulting late-time PDF 
gets modihed as discussed in the Appendix such that the ehective potential gains a quadratic 
term with a well dehned potential, see Eqs. fl55]) and (1551) . While this may be just a curiosity, it 
may contain important physics. Namely, to study properly what happens in the deep infrared 
on accelerating space-times, an RG treatment should be appropriate. There have been some 
attempts to apply RG methods in De Sitter space |l5l4l7l| , but none so far is attempted in more 
general FLRW space-times. While early results are promising, further developments will show 
how important these results are and whether these techniques can be used to obtain results 
that go beyond stochastic inflation, and thus illuminate the curious results obtained in the 
Appendix. 

The late-time solutions we have constructed here can be of help to better understand late¬ 
time dynamics of the inflaton and other interacting scalar helds during inflation, especially 
in the regime when e is not very small. Our effective potential can be used to calculate, for 
example, late-time limit of n—point functions of the scalar held, by making use of Eq. (I28l) . 
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Appendix 


Here we provide a formal procedure for construction of a cutoff-independent late-time PDF 
of an interacting scalar field, whose potential is scale invariant, 


VM = 


(29) 


The cutoff dependence is generated by the fact that free mode functions satisfy an equation for 
which the index u defined in Eq. flTT|) differs from 3/2. To get u = 3/2, one can split the mass 
term in the potential fl2^ as, 

A 

4!‘ 


V (if) = (wi -h 072) ^^^^ , 0 Ji+U 2 = l 


such that z/ = 3/2, or equivalently. 


^2^1^(2-e)(l-6n;iO_9 


is achieved by 


1 3-2e 

uJi = - -e , 


(1-^) 


1 3-2e 

UJ2 = I - - - € . 


6^ 2-e ’ ' 6^ 2-e 

In this case the amplitude of the stochastic force (ITbD becomes cutoff independent, 

(l-e)3 


A ^ Ao = 




(30) 


(31) 


and the late-time limit of Eq. flT^ becomes. 




Poo[^^^ + 6n;2^(2-e)(p - (3-e)e(p 


= 0 . 


(32) 


This equation governs the late-time PDF of the scalar held in stochastic inhation. As in Eq. fj2^ 
above, one can introduce a rescaled effective potential, 

2 


VeS = 


(3 —e)Ao 




(33) 


and the equation of motion for v{(f) = — ln[poo(‘^)/po] is. 


Upon comparing with 


- h" + -t- [U - = 0 . 

, we see that the correct solution is. 


(34) 




( 35 ) 
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From this result we conclude that the cutoff independent solution for the held PDF inherits 
the tree mass term minus a mass term that is proportional to e, such that the effective mass- 
squared is positive (negative) when 2^ —e is positive (negative). We emphasize again that the 
solution fl5S]l is fully cutoff and scale independent. 

While there is no good reason to think that such a modihcation of the theory should correctly 
reproduce early time evolution of the held correlators, the late-time correlators may be correctly 
reproduced. In order to fully understand the implications of the result (05]), one would have to 
perform a suitable RG analysis. 
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